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ABSTRACT 


A [general utility FEH proyrarn v/hich uses isopa:^etric 
elements .and can handle 2- dimensional and 3-diriensional 
stress analysis has oeen developed. The stress analysis of 
spur and bevel gears has been carried cut using this pro- 
gram. Local least squares stress smoothiiv- method has been 
made use of in tlie stress calculations at nodal points. 

Detailed stress analysis has been dene for spur and 
bevel gears. The spur gears have been studied as solid and 
rim gears. The effects of fillet radii are studied for both 
the types. The effects of rim thicknesses are also studied. 
The effects of centrifugal forces on root stresses have 
also been investigated. Detailed stress plots are given 
for each care. 

In case of bevel gears, the effects of fillet radii 
have been studied for tv/o models. The stress distribu- 
tions along the fillet lengths and accross the tooth section 
are presented. For easier application of failure criteria, 
the octahedral shear stress distribution is also given, 

A procedure for bevel g^ear mesh generation is presented. 



CHAPTER 1 


INTRODUCTION 

Gears have found universal application, wherever power 
is to be transmitted without slip between two shafts at con- 
stant velocity ratio. There is no limit to the power that 
can be transmitted. The centre distance is small as compared 
to other types of drives. In the involute tooth form, the 
line of pressure and the line of action coincide on a strai- 
ght line and thus the velocity ratio does not change for 
changes in the centre distance. This advantage has made it 
universally accepted as the standard tooth form. 

1“1- PREVIOUS WORK 

A designer who is set about to design a gear starts his 
design using c^-ntilever beam formula (Lewis -auation) . This 
formula is derived considering gear tooth as a beam of para- 
bolic cross section having uniform strength. Only the load 
acting horizontally is taken into account and the bending 
stress is calculated. The compressive load component is not 
taken into consideration. The compressive stress produced 
by this load reduces the tensile bending stress and increases 
compressive bending stress. The fillet radius plays very 
important role in the stress analysis at the root and these 
are not taken into account in this formula. Jaramillo[l] 


1 
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has made a statement that **i£ the distance of th'e load from 
the clamped e. ^;o is of the order rf mat'nitude of plate thick- 
ness the beam solution ceases tc be valid". Others have also 
commented upon the: inadequacies of cantilever beam theory 
for low bea:.m heiclits. Kov/ever in the absence of alternative: 
approaches, t'ne Lewis equation is still being used for pre- 
diction of stress taking suitable semienperical stress con- 
centration factors. The bending stress alone is not suffi- 
cient to predict failure. In these desicn calculations he 

uses a factor of safety to take care of uncertainty of 

how 

stress calculations. SoA.iph a factor of safety should 
be used? Toe high a factor renders the design uneconomical. 
Another disadvantage is that displacements cannot be calcu- 
lated accurately for different points on the profile. 

Once the design is com^-icte he is faced with questions 
like - 

- Can cantilever bearu theory correctly predict 

the two and three dimensional stress.es actually 
ocBuring in the gears? 

- What are the actual stresses induced and how 

do they compare with the assumed safe stresses? 

- VJ'hat is the stress range at a point'? 

- What are the displacements of various points on 
the profile, so as to take care of noise by 
geometrical ccinpensation. 
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- What is the critical • ecticn and what are the 
stresses th"t cause failure? 

The designer uses one of the following alternatives to 
get answerste the above ouestions. 

(1) Experimental flethods 

[i] Electrical resistance or semi-conductor 
strain gauge measurement 
(ii] Photoelastic method 

(2) Analytical Methods 

(3) Numerical Methods like finite element method. 

Electric strain gauges measure the strain at the chosen 

points. These gauges are difficult to use in the fillet 
regions. Lot of experimental skill and good instruments are 
needed to get reliable results. As this method is not a 
whole field method, results are of limited use only. 

In photoelastic method, first a model is made of trans- 
parent material. The stresses are calculated from isochrn- 
matic and isoclinic fringe patterns. Two dimensional photo- 
elastic analysis is quite simple and straight forward, whereas 
three dimensional photoelastic analysis requires stress 
freezing. This makes the method very involved. Detailed 
post analysis of slices requires experimental skills of 
very high order. 

In analytical methods conformal mapping technique is 
used to represent the gear profile. The selection of 
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transform function needs considerable experience. This ana- 
lysis is quite involved as one has tc try different functions 
to fit the profile. The application of this method t o three 
dimensional analysis loohs almost impossible. 

In the finite element method, the structural idealiza- 
tion of a continuum is accomplished by dividing it into number 
of parts called finite elements. A solution over the finite 
element is assumed initially in the form of a polynomial with- 
out any consideration to the boundary conditions of the struc- 
ture. Using this and the variational statement or differen- 
tial equation of the problem, the typical finite element equa- 
tions (element stiffness, displacement and load matrices) are 
formed. Global equations are formed by adding the elemental 
matrices. Now applying the boundary conditions the displace- 
ments, strains and stresses ar j found at different points. 

This numerical technique is finding larger acceptance because 
of this advantage. These results have proved very reliable 
and are in excellent agreement with the theoretical exact so- 
lutions . 

All the above techniques are being used in gear tooth 
stress analysis. Baronet and Tordion [2] analysed the stan- 
dard tooth form using conformal mapping techniques using 2D 
elasticity and appropriate transform function, the stress 
distribution and geometry factors were found out. 
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V/allace and Siereg [S] have done dynamic analysis using 
FEM. The analysis uses viscoelastic model to describe material 
behaviour for impulsive loading, 

V/ilcox and Coleman [4] applied FEM to symmetric and un- 
symmetric tooth profiles. They compared their results with 
photoelastic results and they were found to be in excellent 
agreement. They also developed a concentration factor that 
takes care of the fillet radius. The stresses calculated 
inside one element were averaged over 10 elements. Statis- 
tical curve fitting was used to extrapolate the stress 
gradient to the surface. 

Chabert, Dangtran and Mathis [5] conducted study on gear 
tooth boundary conditions. They concluded that region bound- 
ed by the lines drawn radially between the adjacent teeth on 
the root circle and l.S-2,eM below the root circle will be 
enough to represent the solid gear for FEM analysis. For 
the worst loading condition when one pair is in contact they 
analysed and plotted the stresses. 

Oda, Nagamura and A cki [6] studied the effects of rim 
thicknesses on stresses. They also concluded that one tooth 
model is sufficient for the analysis. The results were 
verified by experimental means also. 

Aral , Harada and Aida [7] studied the thin rim effect 
on a 9 tooth model. The stress plots are given when diffe- 
rent teeth are loaded. Again the results are compared with 
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experimental results. 

Al”Shareedah and Lehnoff [9] have analysed stralfht 
bevel gear tooth modelling it as a elate. They compared 
these results with tho results of finite element method. 

But their model is of questionable value. The orofile 
of tooth is nowhere near the involute form. The fillets 
have been totally ignored. 

Chang, Huston & Coy [8] have used SAP-IV to analyse 
both solid and thin rim gears. They have studied the effec- 
ts of rim thickness, fillet radii and loading position on 
both these gears. The stress plots are also given. 

1*2. PRESENT WORK 

Most of the previous work done in stress analysis of 
gears by FEM has been on spur gears mostly using triangular 
elements In order to take i„are of the involute profile 
effectively, more nvnber of elements have been used. The 
more number of elements mean more time and pain to prepare 
the data. It is established now that higher order elements 
give same accuracy with fewer elements. It is logical 
to think that isoparametric elements are more suited to take 
care of the gear tooth profile. So it was decided that the 
analysis of tooth should be done using 8 node isoparametric 
element in 2‘-diraensional analysis and 20 node isoparametric 
element in 3-dimensional analysis. 
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Chang, Huston et. ai [8] have studied the effects of 
thin rims for two rim radii. In this study the centrifugal 
forces were ignored. It was decided to investigate this 
aspect and the effects of fillet radii on thin rim gears. 

The bevel gear tooth is subjected to 3o?rfs in all three 
coordinate axes directions, when a single pair is in contact. 
This particular aspect brings doubts in the stress calculated 
by the ccntilever beam (Lewis) formula. The available for- 
mulae are unable to handle this type of loading. Again a 
factor of safety is taken to eliminate the uncertainty. It 
is logical to think that at worst moment of loading, the 
shorter end is more severly stressed. So what is the state 
of stress at this end? What are stress variations along the 
length of the tooth? As in spur gears the fillet radii affect 
the stresses considerably in bevel gears also. So it was 
decided to do three dimensional stress analysis of bevel gears 
using 20 node isoparametric brick element to get a true pic- 


ture of stresses. 


CHAPTER II 


FORMULATION 

Finite element method is used for the analysis of 
all gear models in this work. First the finite element 
method for 3-dimensional analysis is explained briefly. 
Details are available in standard texts [10, 11, 12, 16], 
One with prior knowledge of FEM may skip this portion of 
present chapter. The stress calculations at the node 
points which are not dealt in detail in other texts, are 
given in detail in Section 2-1.3. Mesh generation for 
bevel gears is explained in Section 2-2. 

2-1. FINITE ELEMENT METHOD FOR 3-DIMENSIONAL ANALYSIS 
Potential to be minimized is 

I = / i oedv -/„F. u.dv - / S. u. ds (1) 

Here the first term represents the strain energy, the 
second term represents the work done by body forces and 


the 

third term 

represents 

the 

work done by 

surface forces. 

The 

strain and 

stress 

matrices 

are 
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Strain and stress matrices are related bv [D] matrix as 


{a} = [D] {e} 

Displacements over the finite element are taken 


f ^ 

u 


1n1 

|0| 

lof 


u 

j (ne) 

< V 

> * 

|o| 

1H| 

lo| 


V 

^ E [N] 
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lo| 

|ol 

|N| 1 
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where N are the shape functions. 

Using (5) , strain matrix equation (2) becomes 

{e} = [B] {u}(^®)= [B]i 1^1 

[wJi 

where 


^i , X 

0 

0 

0 

Ni,y 

0 

0 

0 

Ni,2 

^i,y 

Ni,x 

0 

0 

, z 

Ni,y 

Ni,2 

0 

Ni,x 


Substituting the equations (4), (5) and (6) in 
one gets 


(4) 

as 

{u} 


C6) 


equation (1) 


iCe) = 1 ^^,jne ^[B]'r[D][B] dv {u}^° 


/yF[N]'^dv 
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By applying the Rayleigh Ritz method, one gets 


[KJ® 

{u}C’^®) «= {Fg}*^®^ + {Fs}(e) 

(8) 

[1(] 

= / 

V 

[Bl'^[D3[B]dv 

(9) 


= ^v 

F [N]'^dv 

(10) 


_ / 

~ A 

p [N] "^ds 

(11) 


where is the elemental stiffness matrix, 

is the body force matrix and is the surface force 

matrix. Elemental equilibrium Equations (8) are added up to 
obtain global equilibrium equations of the form 

[K] (12) 

where [K] is the global stiffness matrix, is the 

overall displacement matrix and is the global force 

matrix. The displacements are found by solving these simul- 
taneous equations after applying the boundary conditions. 

The stresses and strains are calculated by the use of 
Equations (6) and (4) . 

2 - 1 , 1 Isoparametric elements and numerical integration 

First the finite element with curved boundaries is mapped 
into a cube of side 2 using 


nodes 

2 

= 1 


N • X. 

1 i * 


nodes 

2 Vi. 

i .= 1 


z * 


nodes 

2 


H^Zi 


(13) 


1 


1 == 


i 
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Here shape functions Nj_ may be same or different from those 
used for dependent variables u,v,vf in Equation (5), If 
these shape functions are same as those in Equation (S) , the 
finite element is known as isoparametric element and ha.s 
been used in this work. These shape functions are expressed 
in natural (local) coordinates (r,s,t). Using Equation (13), 
the elemental stiffness matrix Equation (9) becomes 

= Iff dxdydz 

v(e) 

= }/} det[J] dr ds dt (14) 

- 1 - 1-1 

Here [B] matrix is evaluated using the Jacobi matrix [J] , see 
Ref, [10,11,12,16], Jacobi matrix [J] transforms the local 
coordinate derivatives into global coordinate derivatives. 
Using Gauss-Legendre quadrature formula, the above stiffness 
matrix is evaluated as follows 
NGP hJP NGP 

= S 2 Z W.W.W, g(ri,s.,tj^) det [J] (15) 

i=l j=l k=l 

where g = [B]^[D][B] is evaluated at the various gauss points 
and are the weights corresponding to the Gauss 

point (rx,Sj»t]^), These values are ta.ken from the gaussian 
quadrature table [10,11,16], The number of gauss points 
chosen depends on the order of polynomials to be integrated. 
Thus for 3-dimensional brick element with 20 nodes the order 
of polynomial is six and four gauss points are needed in each 
direction for its accurate determination. But it has been 
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found that three gauss point integration is accurate enough 
[10], and has been used. 

Centrifugal Load l^atrix 

For two-dimensional problems, this body force matrix 
becomes 

* ;/ Sl^m [Y]f {p tdxdy (16) 

Here m is the mass per unit volume; is the angular velocity, 

{ is the radius vector from the axis of rotation and t is 
the thickness. Equation (16) after it is mapped into a square 
of side 2 becomes 

* ff [N]'^{y}det [J] tdrds (17) 

4 4 ^ 

This force matrix^evaluated by using Gauss- Legendre quadrature 
formula as explained in Section 2-1.1. In this matrix ip 
is evaluated at the gauss points using Equation (13). 

2-1,3 Stress Evaluation 

The global stiffness matrix is assembled as a profile 
matrix. The system of simultaneous equations is solved by 
.gaussian elimination. From displacements the stresses at any 
point (rj^ ,Si ,ti)^{a} , is calculated by 

ia) = [D] [B(ri,si,ti)] (18) 

Here represents the displacement vector fox that particular 
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element. But it is v;ell established that the stresses 
calculated at nodal points straightaway will be in consi- 
derable error. Experience [10] has shown that the gauss 
points are the best sampling points for stresses. Since 
nodes are most useful output locations for stresses, a need 
was felt for extrapolation of gauss point stresses to the 
nodes in each element by local smoothing [13,1^3. Essen- 
tially in this procedure the nodal values are -Found by an 
exact least square fit to the unsmcothcd values at the 
gauss points. The procedure which is carried out for each 
element is given by [13,14] 

“ngp_ _ NGP__ __ 

i = l i=»l 


NGP _ 

I \Cri,Si,ti)NjCr^,Sj. ,t^) 
i-1 


NGP _ 

i = 1 



NGP _ 

2 Ni(ri,Si,t.3a. 

i= 1 


< : 




NGP 

2 TT C t • s • t ■ 3 ■ 
i= 1 


( 19 ) 
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or [A] = {FS} (20) 

Here the matrix [A] is the smoothing matrix. In the forma- 
tion of this matrix shape functions are one order less 

(i.e. linear functions -f^or quadratic elements) than th^^se 

used for stiffness matrix formation. The right hand side 

OF 

matrix {FS} is the force matrix formed by the use/^unsmoo th ed 
values of stresses at different gauss points. a2» •••> 

are the smoothed nodal stress values at the corner nodes of 
this element. Here HGP is the total number of gauss points 
used in the stiffness matrix formation. This procedure is 
repeated for all the elements. At each node average is taken 
of all stresses given by different elements connected at that 
node. The midside node stresses are calculated by talcing the 
average of stresses at the adjacent nodes. 

2-2. BEVEL GEAR MESH GENERATIOH 

A diagramatic illustration of the bevel gear model is 
Shown in Figure 2-1. For the convenience of construction 
the face of the shorter end is taken as parallel to XZ plane. 
The coordinate -axes are also shown. As shown in Figure 2-2 
the line PQO passes through apex and pitch points of both 
ends. The plane which contains this line is parallel to xy 
plane. Firstly the front view is constructed. Then the side 
and lastly the top views are constructed. 
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Construction of front view ; The terminology is indi- 
cated in Figure 2-3. For a typical involute profile, the 
angle subtended at the centre by the points A and B on the 
pitch circle, (Y), is given by [17] as 

Y = (21) 

Here D is the P.C.D and g is the circular tooth thickness 

Jr 

and is given by 


g 



( 22 ) 


The distance AB on the pitch circle = y x T>^/2 

The angle Yq subtended at the centre by the points C and D 

on the base circle (D|^) is given by [17] 

Yq = Y inv^ (23) 

Here inv is the involute function and is given by [17] 

inv4>*tan<|)-<}) (24) 


where 4).is the pressure angle in radians. 

The distance CD on the base circle ® Yo ^ 

Take any convenient diameter Dj^. If this is the pitch 
circle diameter the pressure angle is given by 



The angle Yj subtended at the centre by the points E and 
F on the circle Dj is given by 


C26) 
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Using Equation (24) invcjjj is calculated. 

Ths distance EF on the circle (Dj) is Yj x D^/2. 

Take few more convenient diameters on the profile and calculate 
the distances on these circles. Measure off these distances 
on both the sides of the axis of symmetry. Complete the pro- 
file by joining all these points. With suitable fillet radius 
draw the curve such that it is just touching root circle and 
the profile. Construct a region bounded by a distance 2M and 
the points on the root circle midway betl^reen the adjacent teeth. 
Using this procedure the profiles of bigger and shorter ends 
are constructed by taking the respective backcone radius as 
pitch circle radius [18]. 

Construction of side and top views : The line PQO is 
taken as the horizontal axis. On this line first the point P 
is fixed. On ths perpendicular dravm at this point mark off 
addendum and dedendum of the bigger end. Point Q is fixed 
after measuring horizontally . . the face width on PQO. 

Measure off addendum and dedendura of shorter end on the perperi- 
dicular drawn at this point. A section )[X’ which is 2M 
distance below the root circles at both ends fixes the lower 
surface of the model. Intermediate lines on which the Z 
coordinates are measured, are drawn depending upon the lengthwise 
sections of the tooth. With the help of front and side views 
construct top view. 
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Construct the FHI'i raoshes on shorter and bigger ends 
taking into consideration the fillet regions* The corres- 
ponding end points of any line on the shorter and. bigger 
ends are transferred from the front view to the other two 
views. The X- coordinates and Y-coordinates of all the points 
lying on this line, are measured from the top view. Simi- 
larly the Z coordinates of the points lying on this line are 
measured from the side view. The uniform load is distributed 
to different nodes on the line RS . Depending upon the angle 
made by line RS in top and side views these loads are 

resolved in all the three coordinate axes directions. 



CHAPTEP III 


RESULTS AND DISCUSSIONS 

As said earlier spur gears (solid and rim) and bevel 
gears have been studied and results are given in this 
chapter. For spur gears, two models SC90 and SC 270 having 
pitch diameters of 90 and 270mm respectively were chosen. 
These gears were modelled as plane stress problems. Effec- 
ts of fillet radii on stress distribution were studied on 
solid and rim gears. Rim gears of different thickness were 
also studied. Effect of speed was studied in thin rim gears. 
Spur gear SC90 was also investigated as 3-dimcns ion al body. 
Next two bevel gears BplOO and BC180 having pitch diameters 
of 100 and 180mm respectively w ere studied. These gears 
were modelled as 3- dimensi onal bodies. Effects o*^ fillet 
radii on stress distribution \fert-; studied. 

For two-dimensional analysis eight node i s opar ame t eri c 
quadrilateral elements and for three-dimensional analysis 
twenty node isoparametric brick elements were used. For 
eight node quadrilateral elements two ^rauss points and for 
twenty node brick elements three gauss points were used in 
each direction for numerical integration. 

Fillet region of the loaded side has been called +ve 
fillet side and the other one as -ve fillet side in the 
present study. 
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3-1. STUDY OF SPUR GEAR, SG-5Q 

This spur gear was studied by [8] and was chosen here 
to check the validity of the prestent computer program)^[jK . It 
may be noted that Ref. [8] used triangular elements with 
three nodes, whereas the present study uses eight noded iso- 

•w 

parametric quadrilateral elements. In the present study gear 
was divided into 28 elements with 111 nodes, see Fig. 3-1. 
Specifications of this gear are as follow: 


Number of teeth, 2 
Module, M 
Ped of gear, Dp 
Addendum, a » IM 
Dedendum, d ® 1.25M 
Load applied, P 

The angle of contact at the top, 

The fillet radii studied, r^ 

3-1.1 Comparison with Earlier Results 

The model S^SO was analysed for tip loading. The rim 
thickness was 0.75M. Figure 3-2 shows the plots of stresses 
obtained from the program. These results are in good agree- 
ment with the stresses plotted by Chang et al [8] for the 


a 18 
= 5mm 
90mm 
= 5 mm 
= 6,25mm 
= 400 kN/m 
= 30* 17’ 

= 0.2H,D.3M,0.4K,0. 


5M 


s ame 
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Principal stress is maximum on the +ve fillet side 

and negligible on the -vo fillet side. Similarly 02 is 
maximum on the -ve fillet side and negligible on the +vo fillet 
side. The compressive principal stress 02 is about 35% more 
than the tensile principal stress aj. Understandably the maxi- 
mum shear stress occurs on the -ve fjllct side. The stress 
gradients are very high near the surface of the tooth. 

3-1.2 Effects of Fillet Radii 

Figure 3-3 shows the stresses across the tooth section 
for different radii 0.2M to O.Sil. Increase in fillet radii 
reduces the fillet stresses as expected. Increase of radii 
from 0.2I'l to 0,5f^ reduces tho principal stresses approximately 
by 35%, The max shear stress oc^curs on the -ve fillet side 
and this reduces approximately by 25% when radii changes from 


0.2M to 0.5M. 

3-2. STUDY OF SPUR GEAR. SG 270 


This gear has PCD Of 270ram. This gear vsis studied as 
solid and rim gear with different rim thicknesses. Effects 
of fillet radii and speed were studied on thin rim gears. 
Specifications of this gear are as follow 


Number of teeth, Z 
Module, M 
PCD of gear. Dp 
Addendum, IM (a) 

Dedendum, 1.25M (d) 

Load applied, P 

The angle of contact at the top, <|)y^ 
Fillet radii studied, r^ 


* 54 
s Smm 
a 270min 
= 5 . 0mm 
a 6,25mra 
» 400 kN/m 
= 25“ !• 

= 0.2M,0.3,0.4M,0.5M 
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As Solid Gear 

Solid gear is modelled as shown in Figure S-lCb) onlv 
one tooth is taken, and a portion of radial thickness of 2M 
below the root circle is included in the tooth for stress 
analysis [5,6], Stress variation across the tooth thickness 
for different, fillet radii is shov/n in Figure 3-4. Princinal 
stresses reduce approximately by 20% when fillet radii changes 
from 0.2M to 0,515. It is observed that the fillet radius has 
little effect on internal stress distribution. 

3“2,2 The Effects of Rim Thicknesses 

Results obtained for this model SG270 arc tabulated in 
Table 1. The rim thickness is varied from ,5M to 3M. The 
predominant principal stresses on +ve and -ve fillet sides are 
taken for comparison. The variation of stresses is negligible 
when rim thickness increases 1 . 5M to 3M. The stresses of 3 . OF'I 
are taken as the base stresses for comparison. When the rim 
thickness reduces to 1 , OM there is a rise of about 15% for O 2 
and beast rim thickness of 0.514 produces almost 80% 

rise in those stresses. The compressive stresses increase much 
more rapidly than the tensile stresses. This observation is 
found to be in agreement with the findings of Chang et\ al fg] , 
3-2.2 Effects of Fillet Radii 

For this model rim thickness v/as taken as 1 , OM and stress 
plot is given in Figure 3-5. 


Stresses at the fillets decrease 
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with increase in fillet radii. The internal stress distri- 
bution is not affected appreciably by the increase in fillet 
radii. The stresses at the fillets reduce approximately by 
25% when fillet radius changes from 0 . 2M to 0.5M. As expect- 
ed the max shear stress occurs on the -ve fillet side. 

3-2.4 The Effects of Centrifugal Forces 

Model SG270 was analysed for 2 different speeds of 600 
and 6000 rpm. The predominant stresses are tabulated in Table 
2, The effects at 600 rpm are negligible. Therefore, the 
results as obtained for the speed of 6000 rpm are only tabu- 
lated. It is seen that even for this sneed, change iri stresses 
is small, of the order of 1 to 3%. Because the centrifugal 
forces are always positive, the compressive fillet stress is 
reduced, while the tensile fillet stress is increased, 

3-3. ANALYSIS OF SPUR GEAR AS THREE DifjENSIONAL ) 30 DY 

Spur gear SG'.)0 was analysed as three dimensional body using 
20 node isoparametric brick elements. Model had 36 elements 
with 286 nodes. Specifications of this gear are the same as those 
given in Section 3-1. 

The predominant principal stresses are plotted in Figure 
3-6Ca) and Figure 3-6(b) for +ve fillet and -ve fillet lengths 
respectively. The 20 value is also plotted on the graph. It 
is interesting to note that the stress value of 3P analysis 
approaches the value obtained by 2D analysis at the middle of 
the tooth. The stresses decrease towards either end of the 


tooth. This is 
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because the points at the two ends are free to move in all the 
directions. The peak values ocour in the middle region of the 
tooth length. Since these peak values obtained by 3-D analysis 
are about 3% more than the 2D values, it is on the safer side 
to analyse the spur gears for uniform load along the length of 
the tooth by 2D analysis. 

3-4. STUDY OF BEVEL GEARS 

Two Bevel gear models BGlOO and BG180 having pitch dia- 
meter's of 100 and 180mm respectively were analysed. The de- 
tails of these models are in the corresponding sections. Each 
model was analysed for four fillet radii. As said earlier the 
stresses of each node were calculated using twenty noded isopara- 
metric elements. The principal and octahedral shear stresses 
across the fillet lengths, tooth sections and tooth profiles 
are plotted in various figures. Shorter and bigger ends are 
identified by SE and BE in the figures for bevel gears. 

Fillet lengths along which the stress variations are plotted, 
were chosen as follow. For any fillet radiuj, the point with 
maximum predominant principal stress was selected on the 
shorter end. The lengthwise distribution is plotted along 
the line joining this point and the corresponding point on the 
bigger end of tooth. 
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The octahedral shear stress distribution is also plotted 
along this line. At each node on the profile of the shorter end 
a perpendicular to the profile is drawn. The corresponding 
principal stresses are plotted on these lines for the distri- 
bution along the profile. 

3-4.1 Bevel Gear, 8G10Q 

This model was analysed for 4 fillet radii of 0. 0.2Mg, 
0,4Mg, 0.55Ms, where Hg is the module at the shorter end. The 
specifications are as given below: 

Number of teeth on pinion 
Number of teeth, on gear Zg 
Module, M 
Pcd of Pinion, Dpp 


Pcd of Gear, Dpg 

Addendum, a = IM 

Dedendum, d = 1.2M 

Face width, b = 6.01‘3 

The analysis is done for the gear 

The module at the shorter end. Mg 

Cone distance, C 

Pitch cone angle, cc 

Addendum angle, 5^ 

Dedendum angle, 6d 
Formative number of teeth Zf 
The load applied, P 


= 16 
= 20 
® Smn 
= 8 0mm 
® 10 Omra 
= 5mm 
= 6mm 
= 30mm 

= 2.657mm 
= 64.031mm 
= 51.34 deg 
» 4.250" (deg) 
= 5.097" (deg) 
= 32 

» 400 kN/m 
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Figure 3.7(a) and Figure 3.7(b) show the variation 
of different stress along the -ve fillet length, and Fig. 
3.8(a3 and Fig. 3.8(b) show the same for the +ve fillet 
length. It is seen that shorter end is more severly stressed. 
For example the predominant principal stresses decrease 
approximately by 35"o at the bigger end than those at the 
shorter end. Also all the stresses decrease with increas- 
ing fillet radii. For example the predominant principal 
stresses decrease approximately by 20 % by increasing the 
fillet radius to O.SS i-Is. Bending stress also shows 

the similar trend. It is observed that maximum of Oy and 
C 3 occurs at about i distance from the shorter end on both 
fillet sides, 

Figure 3-9 shows the variation of stresses accross 
the thickness of the tooth at the shorter end. On the ♦ve 
fillet side is predominant and ^2 is predominant at 

the -ve fillet side as expecti'i. The compressive principal 
stress 02 more than tensile principal stress oj as in 
spur gears. The biding stress shows linear variation 
across the tooth as expected by the cantilever beam theory. 
The fillet radius has no significant influence on the 
internal stress distribution as in spur gears. 
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Figure 3“10(a) shows the variation of octahedral shear 
stress across the tooth section. The maximum value of octa- 
hedral shear stress occurs on the -ve fillet side of the 
tooth and is about 40% higher than that on +ve side. By 
increasing the fillet radius to 0.5H this stress comes down 
by 20% on the -ve fillet side. It is noted that fillet 
radius has negligible effect on the internal stress distri- 
bution. 

Figure 3.10Cb3 shov/s length wise distribution of octa- 
hedral shear stress on the -ve fillet side. The value at the 

bigger end is about 30% less than the value at the shorter 

o-n. 

end and the decrease is mmot^Cv 

Figure 3-11 sho'vs the distribution of predominant princi- 
pal stresses along the profile of the tooth at the shorter end. 
Except for the case ofaerofillet radius, the max values occur 
at points A and S. The values decrease on cither side of 
these points. Thus A-B can be called as the critical section. 
3-4.2 Bevel Gear, EG18Q 

This bevel gear was analysed for 4 fillet radii 
of 0, O.SMg, O. 4 M 5 and O.SMg. Specifications of this 
bevel gear are as follow: 
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Number of teetch on pinion, Zp 

Number of teetch on gear, 1q 

Module, M 

Pcd of pinion, Dpp 

Pcd of gear, Dp^ 

Addendum, & - If* 

Dedendum, d » l.2li 

Face width, b = 6M 

The Analysis is done for the pinion. 

The module at the shorter end, Mg. 

Cone distance, C 

Pitch cone angle, o. 

Addendum angle, 6^ 

dedendum angle, 

Formative number of teetch, 

The load applied, P 


= 18 
» 20 
= 10 mm 
= 120mm 
= 200 mm 
•= 1 0mm 
= 1 2mm 

= 6 OTism 

w 5 , 5 4 mm 
s 134.536mm 
* 41.987 deg 
» 4.250 deg 
= 5.097 deg 
= 24.21 
- 400 kN/m 


Stresses obtained for this bevel gear are shown in 
Figure 3-12 to Figure 3-16. Study of these figures show 
that stress distributions are very similar to those obtained 
the earlier bevel gear, EG190. They differ only in the 
percentage variations. Therefore results are not discussed 

again. 



CHAPTER IV 


COlJCLUSIONS 

Frcra the results obtained for spur and bevel gears 
it is seen that cantilever bear formula (Lewis equation) 
cannot predict the true stresses in the fillet region 
and gives very inaccurate results for thin rim spur poors , 
From the study of 2-D analysis of spur gears, it is veri- 
fied that a distance of 2H below thw- rcot circle is suffi- 
cient to analyse a solid gear. In tv;o spur gears studied, 
it was found that the stresses in the fillet region de- 
creased approximately by 20% when fillet radii increased 
from 0.2I'i to 0.5M. For thin rims, as expected, it was 
found that the stresses increased with decreasing rirn 
thickness. When rim thickness decreased from 2i\ to 0.5M, 
the stresses in the fillet region increased almost by 80%, 
The effects of fillet radii on stress distribution are 
more than those in solid gears. Keeping this in view it 
is recommended that rim thicknesses be at least 0.7SM. 

The centrifugal forces at 6000 rpm, v;hen taken into account, 

did not affect the stresses appreciably. Hence it is 

sty ess 

justified to neglect these effects for the rcot ^cal cul a- 
tions . 
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Three dimensional analysis of spur gears shovrsd 
that the stresses in the mid region of tooth reach the 
2~ dimensional stress values. At both ends of tooth 
the stresses are lower than those at the midregicn. 

Hence the 2-dimensi cnal analysis of uniformly loaded spur 
gears is justified. 

In case of bevel gears it v/as observed that the 
shorter end is much more scverly stressed than the bigger 
end. In the two bevel gears studied, it was observed that 
stresses in the fillet region decreased by about 20% at 
the shorter end when fillet radii changed from 0.0 to O.SMs, 
This decrease was less at the bigger end. The stress 
distribntion accross the thickness v;ere as in spur gears. 

The octahedral shear stress was also maximum at the 


shorter end. 
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Fig.3'10 (o) Octahedral stiear stressi#it across 
the tooth section, BG 100 
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